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Abstract 

Brownian motion is modelled by a harmonic oscillator (Brownian particle) 
interacting with a continuous set of uncoupled harmonic oscillators. The 
interaction is linear in the coordinates and the momenta. The model has an 
analytical solution that is used to study the time evolution of the reduced 
density operator. It is derived in a closed form, in the one-particle sector of 
the model. The irreversible behavior of the Brownian particle is described by 
a reduced density matrix. 
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I. INTRODUCTION 



In previous works |]-|J we have studied the time evolution of a quantum oscillator cou- 
pled to a dense, but discrete (finite) bath of harmonic oscillators. For such system an 
irreversible behavior has appeared as a consequence of averaging in time the evolution of 
the characteristic quantum-oscillator variables (macroscopic quantities), since time evolu- 
tion splits in very different scales: One related to small fluctuations, which are erased by 
averaging, another one related with recurrence phenomena, which are far enough of labo- 
ratory observational times, and the last one connected with observable phenomena, which 
involves irreversibility. 

In Ref. [0J it has been shown how to pass to the continuous bath. A resonance, which can 
be isolated and leads to an evolution for a macroscopic period of time (the same period as 
in the discrete case), has arisen because of the energy of the quantum oscillator is embedded 
in the continuum. One particularity of this continuous limit is that two of the three time 
scales have become irrelevant since they have zero measure with respect to the remaining 
time scale. However, this last scale is not exclusively governed by the resonance (associated 
with an exponential decay) but also by contributions coming from the semibounded feature 
of the energy spectrum. 

Leaving these contributions aside the main behavior of the quantum oscillator is an ex- 
ponential decay towards the equilibrium with the bath, and can be described only with 
the contribution due to the resonance. The present work is the conclusion of the previous 
papers |5|-f|, where a novel method was used to work directly in the continuum, including 
the exponential decay law in quantum mechanics. In this paper we continue with the devel- 
opment of the formalism we have called "Minimal irreversible quantum mechanics," where 
time asymmetry can be represented through the choice of a subspace of "admissible" or 
"regular" solutions of the evolution equation. 

The main idea goes as follows. Usually rigorous quantum mechanics must be formulated 
in a Gel'fand triplet Q 
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ScHcS x , (1) 

where: 

•S is the space of "regular states" or test-functions space, corresponding to Schwarz-class 
wave functions, that are considered as the "physical" states. 

•H is the space of "states," or Hilbert space, introduced to extend the notion of proba- 
bility to a larger space and to use the well-known spectral theory of Hilbert spaces. These 
states correspond to square-integrable wave functions. 

•S x is the space of "generalized states", or rigged Hilbert space, namely the space of 
linear (or antilinear) functionals over S, which are essentially used to find the spectral 
expansion of the regular states (e.g. Fourier expansions). 

Let K be the Wigner or time- inversion operator. As usual the evolution Hamiltonian H 
is time symmetric, i.e. 

KHK^ = H. (2) 

In the wave function representation the action of K coincides with the complex conju- 
gation, so it is defined over S by 

K<p(x) = <p*(x). (3) 

Thus 

K : S -> S. (4) 

Therefore S is also time symmetric. 

But the real universe and macroscopic objects have clearly time-asymmetric evolutions. 
Therefore the task of this paper and the preceding ones is to explain how this time asymme- 
try appears while the quantum mechanical laws of the universe (embodied in Ti) are time 
symmetric. The usual and successful explanation is based on coarse-graining: Macroscopic 
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objects have a huge number of dynamical variables and we can measure and control only 
a small number of them, the so-called relevant variables. If we neglect the rest of the vari- 
ables, the irrelevant ones, we obtain time-asymmetric evolution equations. Nevertheless in 
paper || (according to the line of thought pioneered in Refs. 0, |]I01 , fTTf ) we have follow a 



different way. We have developed a sort of minimal irreversible quantum description, which 
reproduces time asymmetry from the basic microscopic level directly, where the key point 
is the presence of resonances (and additional hypotheses we have extracted from Ref. |§). 

Obviously we want to obtain the standard results making minimal changes to the well 
established and usual quantum mechanics. If we change Eqs. (0) or (§), it is almost sure 
to find experimental problems. So the minimal modification is to change Eq. (^) defining a 
new test-functions space $ + C S such that 

K : $ + -»■ ^ $+. (5) 

In this way K is not even defined over the space of regular states $+ and a time-asymmetric 
evolution arises. 

This can be done if we postulate, as we have done in Ref. ||,[] that all the "regular" or 
"admissible" states belong to a space Ti + ~ #(^+) and also to S. Then $ + ~ 9(H^_nS) [the 
time inverted states belong to a space 7i- ~ #(if 2 ) and $_ ~ #(if 2 HiS), respectively] , where 
9 is the Heaviside step function that gives the restriction to the positive real energy axis and 
H± are the Hardy class function spaces |J .0 

An "irreversible" quantum theory based on a Gel 'f and triplet 

$± C H± C (6) 

is feasible and it yields physical results, as the dominant experimental decay of unstable 
states, if the test-function space $ + is so chosen. This will be valid for systems where the 



x This postulate has been motivated in cosmological-global considerations in Refs. and 
2 As spaces Tt- and 7i + are isomorphic they are normally called TL [14]. 
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existence of resonances dominates the evolution for the relevant period of observational time. 
We have shown that, what it is done in the quoted papers [|9], PH|| , and is essentially a 
minimal modification of the ordinary reversible quantum theory. In fact, from now on we 
will consider that: 

•$ + is the space of "regular states" or test-functions space, that are considered as the 
"admissible" states. 

•TL+ is the space of "states," or Hilbert space. These states are again particular square- 
integrable wave functions. But in paper || and in this work we consider that only $ + 
contains the "admissible" states. 

is the space of "generalized states," or rigged Hilbert space, namely the space of 
linear (or antilinear) functionals over $ + , which are essentially used to find the spectral 
expansion of the "regular states" in Sec. III. 

The spaces with subscripts "— " contain the time-inverted states of the corresponding 
spaces with subscripts "+" . 



Friedrichs model |2S] was studied using this approach. In this work, we show that this 
idea can be used to take a slightly different point of view in studying dissipation phenomena 
of quantum Brownian motion. This more complex model will force us to generalize the 
definition of space $+ although the roles played by the characters in the triplets $± C 
7Y± C will remain the same. 

Brownian motion has been extensively studied in the literature (we will only quote those 
papers particularly relevant to our line of work). E.g., in Ref. jHJ it was shown that for a 



system composed by a finite number of linear interacting oscillators a dissipative behavior 
can be found in the limit of a dense system (continuous spectrum). But, in this work 
we are concerned directly with dense systems with continuous spectrum. The presence of 
this continuous spectrum allows us to study the decay processes using analytical properties 
familiar in scattering theory H|T(J . The model is a very well-known and widely used system, 
consisting of a harmonic oscillator coupled to an infinite and continuous bath. In this paper, 



as in Refs. [15,17|, the bath is composed by an infinite collection of harmonic oscillators 



and the interaction is modelled to be linear and characterized by the spectral weight, but 
otherwise arbitrary. We show that the oscillator reaches a final equilibrium state via a 
damped evolution which is mostly exponential. We also show that some deviations from 
this exponential decay law (for very short and very long times) appear, which are intimately 
related with the presence of a lower bound of the energy. 

In Sec. II the whole system (single oscillator plus the bath) is described and the Hamil- 
tonian is introduced. 

In Sec. Ill we diagonalize (in normal modes) the Hamiltonian. In the process of diago- 
nalization some problems emerge, such as the lost of the discrete part of the energy spectrum 
|28| , |T8|1 . We can bypass these problems, if we use our definition of "regular" states. Then we 
can perform an analytical continuation of the spectral decomposition of the Hamiltonian, 
promoting the energy to complex values. To reach a successful interpretation of the results 
we require to generalize the definition given in paper || to the model we are now studying. 
The mathematical bases of this generalization are shown in Appendix B.Q This appendix 
also contains a rigorous mathematical understanding of the problem. 

In Sec. IV mixed states and their evolution law are considered. 

In Sec. V we deal with a very particular initial condition: An oscillator in a zero- 
temperature bath. We find the reduced density operator and show that the equilibrium 
state is reached. We accurately describe the time evolution of the system and estimate the 
Zeno [19| and Khalfin |2(J effects for very short and long times, respectively. These are the 



deviations from an exact exponential decay law. Finally we show that our solution satisfy a 
Lindblad master equation when discarding these deviations from the exponential behavior. 

Finally in Sec. VI some questions concerning irreversibility, already considered in papers 
2l| and ||, are discussed. 



3 The reader who is not familiar with rigged Hubert spaces and functional analysis can see Refs. 
|] and @. 



6 



We state our conclusions in Sec. VII. 

Three mathematical appendices complete this work. 



II. PARTICLE-BATH MODEL 



The system is a Brownian particle represented by a harmonic oscillator with natural 
frequency Q. It is well known that for a finite bath it is not possible to prove convergence in 
an equilibrium state in the limit t — > oo because of the existence of recurrences [|I^,E3-p6 



However, for large systems these recurrence times become extremely huge and we can elim- 
inate them by passing to the limit of an infinite continuous bath. Therefore, in this paper, 
we consider the oscillator in contact with a bath, already modeled by a continuous set of 
harmonic oscillators with natural frequencies to. The coupling between the system and the 
bath is assumed to be linear with strength g(u). The Hamiltonian for the composite system, 
in terms of creation and annihilation operators, is 

roo poo . v 

H = Via) a + J duub]J) w + X J duog(uj) (a^b^ + b^aj . (7) 

The first term corresponds to the system, the second to the bath, and the third one corre- 
sponds to the interaction between them. In order that the Hamiltonian would be positive 
definite we must require pT|J27j1 that g(0) = and 

Q>A 2 fdo/M. (8) 

JO UJ 

This is an important condition which selects the kind of spectral densities appropriated to 
lead to an irreversible evolution. For example, the ohmic case which is frequently used in the 
literature must be disregarded, unless a cutoff is used. (Operators b w and b^ are rigorously 
defined in Appendix A). 

The Fock basis is the tensor product of the Fock basis of the isolated harmonic oscillator 
and those of the bath, namely 



\n,oj-i...u m ) = \n) ® lut-.-Um) , (9) 

where |u;i...a; m } represents a state with m quanta in the bath, each one with frequency Uj 
U = 1, ...,m). 

The total number of quanta is conserved allowing us to solve the problem by sectors 
(block diagonalization) . The one-particle sector is referred as Friedrichs' model |28| and 
contains the relevant information that we need to compute physical quantities (see Sec. V). 



III. NORMAL MODES OF THE HAMILTONIAN AND ANALYTIC 

CONTINUATION 



The linearity in the coupling term of H allows us to easily find a new set of uncoupled 
harmonic oscillators (normal modes), such that 



/ = / du b%, (io) 

Jo 

H = / duubt,b u , (11) 







where 

b u = Ca+ I du'Q^bu'. (12) 



From a straightforward calculation |6, HJ21], using the Heisenberg equations of motion, we 



obtain the coefficients of the unitary change of variables which diagonalize the Hamiltonian, 
precisely 

*„(u')=5(u-u')+ *» 9 V\ (13) 

[U — UJ 1 + IE) 

and 
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where 

oi(z) = z-n-X 2 [°° du^-. (15) 

JO Z — OJ 

This function, which is the inverse of the reduced resolvent of H in the one-particle sector, is 
not entire because it has a cut along the positive real axis corresponding to the continuous 
spectrum of the Hamiltonian. If a(z) 7^ for all z G C, except for a possible real and 
negative u>q such that oe(u)o) = 0, an isolated solution appears, which is non-analytic in 
A . We do not consider this case henceforth, since we are interested in analytic solutions 
satisfying condition (||). 

If a(z) = has no real solution it is not possible to find an operator a, such that a — > a for 
A —>■ 0. In this case we have lost the particle number operator corresponding to the discrete 



part of the spectrum of H and we do not have the correct form of H when A — ► |30| . This 
problem can be solved promoting the energy (or frequency) u to be a complex variable z. 
We define (3(z) = [a(z)]~ . It can be proved that f3(z) has the same analytic structure than 
the one of the coefficient S(z) of the scattering matrix ||. f3(z) is a meromorphic function 
on a double Riemann sheet with a cut along [0, +00). /3±(u) = [a±(u>)]~ 1 = /3(u ± is) are 
defined on the upper and lower half-planes of the first Riemann sheet R[ (physical sheet), 
and have meromorphic continuations to the lower and upper half-planes, respectively, in 
the second sheet Rn (unphysical sheet). For simplicity we consider g(z) such that the 
analytic extension of /3 + (z) into the second sheet has a simple pole zq = ojq — 1 7 [7 > and 
a + (zo) = 0] in the lower half-plane. Also (3-(z) has a simple pole z$ on the upper plane in 

Rn- 

We can now study the meaning of Zq. From the role played by zq in the evolution equation 
we know that (Im^o) -1 = 7 _1 is the mean life time of the unstable state = \0,v) 

and (Re^o) is the shift of the bare frequency Q [see 0] and also Eq. ( fKp l. But, z is the 
root of a+{z) and from Eq. (|T5| ) we can estimate, up to the second order in A, 
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z = n + \ 2 P [ duj^-^- - m\ 2 g 2 (n), (16) 
o 

where P denotes the Cauchy principal part of the integral.^ The mean life of the unstable 
state and the shift frequency are given by 

7 = 2vrAV (fi) (17) 

and 

5n = \ 2 P [dcu^-. (18) 
J i ' — UJ 

o 



which are well-known expressions in the theory of unstable systems ||31|| , usually derived 
from the Fermi golden rule. 

Regarding the coupling function of the form g(u>) ~ uj n we find that the ohmic case 
(n = 1) without cutoff does not satisfy the positivity condition (|j). If we call 7i/ 2 the 
coefficient for the subohmic case (0 < n < 1) and 72 the coefficient for supraohmic case 
(n > 1), it is easy to prove that 

72 < 71 < 7i/2- (19) 
Now we will find a generalized partition of the identity / and a generalized spectral 



decomposition of H that recovers the discrete part of the spectrum |10|j21|j29|] . In order to 
do this let T be the curve of Fig. 1. It lays on Ri for (3-{z) and on Rjj for f3+(z). We define 
the analytic function of z G C 

a r (z) = z-n-\ 2 ( dz'^\, (20) 
Jr z — z' 



^The principal part comes from the well known identity between distributions 

1 1 

P iir5(x), x e R 



x + IE x 
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which generalizes Eq. (|15|). 

To find the partition of the identity and a expansion of H we will use in some adequate 
analyticity properties.^ Thus we will define a space $ + C H + of states \ip) such that 
the function (0, uix, LO n \<p) = <po(uJi, uj n ) would have an analytic continuation, for each 
variable u { (1 < i < n) to a region that include the singularity z . This space <£> + would 
be our space of "regular," "admissible" or "physical" states. Precisely, generalizing what 
we have done in paper f5j, we will chose $+ such that its states would satisfy Eq. ( |B3"D 
of Appendix B. Analogously, the space $_ of the "unphysical" time-inverted states would 
satisfy Eq. ([B4|) of that appendix. With this choice the analytic continuation that we will 
perform has a rigorous meaning, since the operators act in a space $ + which endowed with 
adequate analytic properties. The demonstration of this fact is a mathematical problem, 
which is considered in Appendices B and C, where we generalize previous results. Then if 
\ip) G $+ and G from Eqs. (|T0D and (|TT| ) using appendix B and following the similar 
demonstration of paper 0, it can be proved that 

= M(5<->*5<-> + f dz b<-WP)\<p), (21) 



H\<p) = (^{(zoa^a^ + dz zb^H^)\<p). (22) 



r 



The residue at zq contributes to the first terms of the r.h.s. of these generalized partition 
of the identity and spectral decomposition of H, as in paper [10[|, and, in a weak sense, the 
two previous equations can be written as 

/ = s(->*a(-> + [ dz bi~H^ (23) 



r 



H = z Q a^>a^ + J r dz zb^Vf*, (24) 



'Properties of this kind were already introduced in previous works ^J^,|TT|,^] . 
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The creation and annihilation operators in all these equations reads 



r(-) 



a + A dz- — o 



(25) 



-6t 



and 



The distribution 



bt ] = b u + 



bt» = hi + 



means 



\g(u) 



a_(a> 



fl t + A r dz^^ 

JO [Zq ~ Z\ 



a + M dJ MLl 

o to — u' + ie 



i f + A / dJ 
Jo 



9W 



0J — OJ — IE 



b,„> 



(26) 



[*0 -z] 



-duj 



-dz 



-da; + 2mf(zo) 



(27) 



/o [zo — w ] + Jr -z — -2 -'o Zq — 
for every well-behaved analytical function f{z). Observe that b^* does not change if we 
replace / r by / °° because, in this case, no pole is crossed, since ce-(z) has not poles (in Si). 
Nevertheless b)~' does change because the pole is crossed when we modify the integration 
contour. We have shown this fact explicitly putting rj + (uj) in place of a+(u), where 

1 



(28) 



rj + {u) a+{u) a' + (zo) 
being 5(z — z ) the extension of the Dirac delta defined d la Gel'fand and Shilov As a 

consequence of these facts a^* ^ d^~^ and b^~>* ^ ■ The star operation corresponds to 
the analytic generalization of the complex conjugation, which, acting on an analytic function 
f(z), is defined by| 



(29) 



6 It corresponds to the symbol $ of paper ||. Here we follow the notation * of paper [21] in which 
we have also studied this model. 
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From Eqs. (|25|) and (|2q) we see that we have four annihilation operators due to the 
presence of complex eigenvalues of H with the corresponding doubling of solutions, since 
we have a pair of complex conjugate values. They are generalized eigenvalues of the two 
analytic continuations of H into the lower (— ) [upper (+)] -complex plane. These operators 
are 

S<->; «W; &H; &« 

The vacuum is the state annihilated by any annihilation operator. The Bogolubov transfor- 
mation of Eqs. fl25|) and (p6|) does not mix creation and annihilation operators, therefore the 
vacuum just defined is actually the same state defined as the vacuum of the noninteracting 
system+bath. So from Eq. (§) the vacuum is the state |0) <8> \v) = \0,v), where \v) is the 
vacuum of the bath. 

The corresponding creation operators are 

aH* ; 3 <+K &H* ; 

Starting from the common vacuum, by applying successively the operators a*- - -** and 
b^*, the Fock basis {|~)} is built, and with 5^* and we build up the Fock basis 
{| )}. In the case of a' - )* and a^* the corresponding vectors in the Fock bases of the 
one-particle sector are generalized eigenvectors of H with purely complex eigenvalues. They 
represent unstable states, i.e. a^*\0,v) = l,v) is a one-particle generalized eigenvector 



of H corresponding to a complex eigenvalue z and a^* 0,vj = l,v) is a one-particle 
generalized eigenvector of H corresponding to a complex eigenvalue z^. In this way we are 
able to develop a second quantized version of the theory of unstable states |^T] . 



Now we have two different number of quanta operators, = J °° dtob^*b^ and 

iV( + ) = / °° dub , which are not Hermitian. So two different Fock bases can be built 
satisfying 

iV ( ~ } |n, uJ^7..u m ) = m \n, uhT..u m ) , 

(30) 
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\n,u>i...uj m ) — m\n,u>i...uj m ) . 
The spectral decomposition of the Hamiltonian reads 



00 00 poo POO 

H^> = V V / ... / dux...du m (zqu + lox + ... +u m ) 

x |n, uhT..u m ) (n, Ui...u m \ , (31) 



n=0 m=0 ' 



which acts on the right of the Fock space generated by basis {|~)}. But the "same" Hamil- 
tonian can also be written in the following way (using the other analytical continuation, in 
which case it is evident that the next equation is only weak, and it has a precise meaning 
operating between \<p) e <E>_ and G $+) 

°° 00 poo poo 

= V V / ... / dujx...duj m (zqU + Ui + ... + u m ) 

trz n ±r n Jo Jo 

x \n,ui...(jj m } (n, u^T..u} m \ , (32) 



n=0 m=0 



which acts on the right of the Fock space generated by basis {| )}. In the same way, the 
identity reads 



^ ^ poo roo 

= V V / ... / du 1 ...du m \n, uf^.Um) (n, Ux...u r , 

Z~ n -T^n JO JO 



n=0 m=0 



(33) 



^ ^ /'OO /"OO 

= y2J2 ■■■ dui...du m \n,ui...u m ) (n,cJ^7.u m \ . 



n=0 m=0 ' 

As the eigenvalues of Eqs. (|31|) and (|3^) are complex, in order to deal with unstable 
states we must find an adequate mathematical structure beyond the Hilbert space. In 
fact, these states are generalized states. Thus, in the Appendix B we see that kets |~) 
and | ) are well defined in a rigged Hilbert space formalism, i.e. they must be thought 
as antilinear functionals acting on test spaces <3?± and, as elements of a vector space, they 
belong to the duals of <3>±, symbolized by $±. They define a double Gel'fand triplet structure 



$± c H± c IH^. 

What about energy conservation? The trouble emerges because the eigenvalues of the 
Hamiltonian are now complex; thus some states decay in time (e.g. vectors of which 
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vanish for long times, see Sec. V). Since some states vanish we may ask ourselves how the 
conservation of energy can be possible. The answer is that energy is conserved anyhow. In 
order to demonstrate this fact we will calculate the mean value of the Hamiltonian in a state 



oo oo 



\<p(t)) = E / — / dto 1 ...du m c n {(*)x...LO m ) \n,uJi...u m ), precisely 

nm o o 



E=(<p{t)\HU\<p(t) t 

oo oo 

= X! H /•••/ duj 1 ..du m du[..Ju f m ,cl(u 1 ...u m )c n >(u[...uj f m ,) 

nn' mm' q q 

xe -i{z n'+^' 1 +...+^' m ,)t e i(z^ 1 +...+u Jm )t^ Qn + Ui + + ^ 

x (n,ui^..uj m \n',uj' 1 ...uj' m ,) . (34) 



Taking into account the orthogonality relations ( |B7| ) and ( |B9| ) shown in Appendix B, 
and that z = u — z^, 7 > 0, Eq. fl3H| ) reduces to 



00 00 



E = J2 ■■■ du x ...duj m \c (uJi...uJm)\ 2 (wi + •■■ + w m ), (35) 
m 

which is time independent. Thus energy is conserved. Conservation of the norm and the 
number of particles can also be demonstrated in an analogous way changing H by /. 

Finally we can observe that Eqs. (|B7| ) and (|B9|) show that the generalized eigenvectors 
have null norm and energy (with the exception of those with n = 0) |5|,|7|, |2l| , |34]| . In the 
literature they are called Gamow vectors, they are generalized states, and they represent 
just idealized mathematical states (see Appendix B), as is the case of the plane waves. 



IV. MIXED STATES: ITS EVOLUTION 



A general pure state belonging to $ + (see Appendix B) can be written as 

/•OO /"OO , . 

|*) = VV/ ... / (kj!...(kj m c n (u 1 ,...,u m )\n,u 1 ...u m ), (36) 
- - Jo Jo 
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where |n, Ui...u n ) G and the most general density^ matrix can be written as 

roo roc 

p = VV/ ... / dux...du m du[...du' m ,c nn i (ux,...,u m ,u[,...,u' m ,) 
, , Jo Jo 

nn' mm' 



x\n,u l ...u m )(n,uj' 1 ,...,uj' ml \. (37) 

If this p is the initial state p = p(0), the evolution law of pit) reads 

p{t) = e- iH{ ' h p(0)e iHi+h . (38) 

As H is only self-adjoint in a generalized way acts in a different way than = 
H^~^ and there are right and left eigenvalues, 

H^\n,ui...uj m ) = (z n + ui + ... + u m )\n,u 1 ...u) m ), 



(39) 



{n,ui...u m 

Then we have 



= (zn-n + uj\ + ... + u m ) (n, uj\...uj ri 



p / t ) =y e -¥n+n')t e -iMn-n')ty du 1 ...duJ m dw'..dw', 

, , Jo Jo 

nn' mm' 

xe -(-i+---+- m )* e iK+---+^)* Cnn ,( Wl; ... )a ; m ^; ) ... )a; ^^ 

x\n,Ui...u m )(n',u' 1 ...u' m ,\. (40) 

For an arbitrary initial state p(0) a time dependent asymptotic (£ — > +oo) state is 
reached. The explanation of this fact is simple. The modes of the bath are independent of 
each other [see Eq. (0)], and so we cannot expect that the bath reaches equilibrium (cf. Ref. 
HI). Thus 



7 More accurately, we would say, that it is the most general possible decaying density matrix, as 
we will see. 

8 Recall that H is self-adjoint in the Hilbert space, where H = 7i x , but in the generalized Hilbert 
space this property essentially becomes Eq. (|39|). 
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rco poo 

Pit) -»• P*(t) = E / - / d0J 1 ...duJ m duj' 1 ...duj' fe -i(^+-+^m)t^ 1+ ...+u> m ,)t 

xc 00 (o;i, ...,w m ,w[, ...,w' m ,)\Q,wi...w m ) (o,w , 1 ...u/ m , . (41) 
For completeness we also write the evolution equation for the density operator, 



^ -_^ e -i(n+n')* e -io, o(n-n')*^ / / d U 1 ...du m du[...du' m , 



co poo 



dt nn> mm> J ° J ° 

xe -i(u Jl+ ... + u Jm )t e i { u J [ + ... + u J ' ml )t^ + Ui+ +Um _ z * n , _ w / _ _ ^ 

xc nn '(wi, - , w m , wi, w m ') |n, W!...o; m > (n', a;}...w^ , (42) 
which is clearly equal to 

^ = ~i (H^p - pHW) = -iLp, (43) 



where L is the generalized Liouvillian operator ||21|| . So we see that the density operator 
follows an evolution described by a generalized Liouville-von Neumann equation. 

In spite of the result obtained in Eq. ([IT]) in Sec. V we show that the reduced density 
operator p r , which is obtained by taking the partial trace with respect to the environment 
modes, reaches equilibrium, namely a time independent state. An equivalent way to find 
the equilibrium state, closer to the spirit of our formalism, is to use a particular space of 
observables, as in paper 



V. REDUCED DENSITY OPERATOR 



As an illustration of the formalism we consider a simple example where the initial state 
is a very particular state of the composed system, 

p(0)=ps(0)®p E (0), (44) 

where 
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p s (0) = cu |1) (1| + cio |1) (0| + coi |0) (1| + coo |0) (0| 



(45) 



is the initial state of the discrete oscillator, or the initial reduced density operator, (with 
en, c 00 > 0, c u + c 00 = 1 and c w = c^) and 



MO) = \v) (v\ 



(46) 



is the initial state of the bath, which does not have any quantum, namely, it is in the 
ground state. This corresponds to a bath at zero temperature T = (in the zero and 
one-particle sector). The main features at any T can be reproduced but we begin with this 
example because the mathematical computations are easier (recall that this model can be 
decomposed in sectors of constant number of quanta). Also our initial conditions are such 
that there is no correlation between the oscillator and the bath. 

Our aim is to find the time dependence of the reduced matrix elements. It is derived 
from the time evolution of the density operator 



p{t) = e- lH{ - h p(0)e lH{+h = e - iffHt /(-)p(0)/W e ,ff(+)t , 
where 1^ and 1^ are the identities in spaces $ + and respectively^] [see Eqs. 

P(0) 



(47) 



P(t) 



-iZQt 



l,v)(l,v 



-iuit 



x I e l2 °* 



1,v)<1,t; 



/■oo 

+ / duo e 
Jo 

poo 

+ / duo 1 e w 1 
Jo 



0,u}) (0,u 



0,o/> (0,W 



(48) 



We have only considered the terms of the identity that correspond to zero-particle and 
one-particle subspaces, since, from the conservation of the number of quanta, there is no 
contribution of other terms. We emphasize that no approximations were carried out up to 
here. 



9 The difference in the conventions with respect to the use of + and — is the following. In operators 
+ and — are related with the analytic continuations for dbie, while + and — in spaces are associated 
with the time evolution, which is only well defined for positive or negative times, respectively. 
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Once we have the time evolution of the density operator, the following step is to get the 
reduced density operator, tracing over the basis corresponding to the environment, 

00 roo rod 

pJt) = tr E p(t) = V / ... / dui...du m du' l ...du' m 

m=0 

x (u 1 ...u m \p(t) \u[...u' m ) S(ui - Lu[)...5(uj m - u' m ), (49) 

where the m = term means (v\ p(t) \v) . 

As we have said, the contribution of terms m = 2,3, ... vanishes. Therefore, using the 
relations between "new" and "old" bases [Eqs. (E3]) and (Bfi|)] and the conservation of trace, 
we obtain a positive definite reduced density operator: 

Pr (t) = C U P(t) |1> (1| + C 10 A (t) |1) (0| + 



c iA*(t) |0) (1| + {c 00 + c u [1 - P(t)}} |0) (0| , (50) 

where 



A o(t) = + / du e-** ,\ K (51) 

and P(t) = |A (t)| 2 is the survival probability of the state with only one quantum in the 
discrete part. 

We can write P(t) as the sum of four terms where the first one, j^f^p , shows an exact 
exponential behavior. Expanding |a^(,2o)|~ 2 as 1 + 0(A 2 ), we split the probability into 
two terms, one containing the purely exponential contribution and the other that we call 
"background," giving rise to derivations from that purely exponential decay law, so that 

P(t) = e~ 7i + background. (52) 

If we take a time neither very short nor very long, the background will be smaller than the 
purely exponential term (for A 1) and it can be neglected, which leads to an exponential 
decay-law. This is not true for short times since tt(0) = 0, which leads to the so-called Zeno 



effect Ui~0 |. For very long times the exponential term will decay faster than the background 
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will do, which is known as Khalfin effect ||20|| . We can force P(t) to have an exponential 
appearance by defining the decay rate T(t) to be time dependent, namely 

P(t) = e~ m \ (53) 

with 

F(t) = 7 In (l + e 7 * background") . 

Obviously for an intermediate time the background can be neglected and T(t) ~ 7. The 
main restriction, imposed by the Zeno period, is 

^(0) = -r(0) = 0. (54) 

For a very long time the decay probability has also a non-exponential contribution as a 
consequence of the semi-finiteness of the energy spectrum. From papers |||7| we know that 
the survival amplitude goes to zero as t goes to infinity as a consequence of the Riemann- 
Lebesgue theorem. Then the behavior of A depends on the small-frequency behavior of 
g 2 (oj). For small cu, ct+(u>) ~ — u> , where condition (|j) was also used. Then the form 
of Ao depends essentially on g(u) for large t. As an example we consider the case where 
g 2 (co) ~ u; n exp (— ^-J , where A is a cutoff [see paragraph before Eq. (|19|)1. By evaluating 
the survival amplitude we have 

A (t) = A 2 [~ du-^pLe-** ~ f^dco W »e-*-*exp (~) , 
Jo Jo \ A / 

where the contribution of high-frequency terms is negligible. Performing the change of 

variables out = x, we obtain 

A (t) ~ jf 1 dx aTe* (l - ^ + ^ + . (55) 

We can see that the survival amplitude merges into an algebraic long-time tail. The first 
relevant contribution behaves as t~ n_1 . As a consequence the decay rate for long times must 
behave as 
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T(t) oc ^. (56) 
The behavior at short times and intermediate times coincides with those obtained in Ref. 



1711 . 111 Fig- 2 we show the qualitative behavior of P(t). Zeno's time, tz, and Khalfin's time, 
tx, are not in scale in the picture in order to show the three different contributions to the 
decay probability. 

Eq. fl5"0|) is the exact solution to the proposed problem, without taking any approxi- 
mation. The first, second, and third terms will vanish for t — > oo; in fact P(t — > oo) = 
and the same happens for A (t) (recall that P = |A | 2 ). The first term of ( |5"ID has the 
factor e~% t and the second one will tend to zero because of the Riemann-Lebesgue theorem. 
The probability of having the vacuum will grow in time. This means that all quanta in 
the discrete spectrum, except the ground state, decay into the continuum. So we find the 
equilibrium reduced density operator 

P* = (coo + c u ) |0> (0| = trp |0> (0| = |0> (0| . (57) 

As expected, the equilibrium state is the vacuum, namely for t — > oo there are no quanta in 
the discrete spectrum, because the initial quantum has decayed into the bath (the discrete 
oscillator has spread its energy over the infinite oscillators of the bath with a distribution 
centered at the shifted frequency uq) P,|T5|,|6| . This means that the discrete harmonic oscil- 



lator has thermalized at T = 0. A similar result was recently obtained in Ref. 

In order to check the compatibility of the solution found in Eq. fl5"U| ) we first briefly sketch 
the main points of the derivation we have done. We have obtained the exact solution of the 
Liouville equation. As a particular case, we have considered an initial condition restricted 
to the zero- and one-particle sectors and we have traced this solution over the environment 
modes. In that case, the survival probability P(t) in Eq. (|50l) can be approximated by 
an exponential behavior, when the background contribution is neglected. Now the solution 
for p r {t) obtained through this approximation can be derived from a master equation of a 
Lindblad's form, where the Lindblad generator is proportional to the destruction operator a, 
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since in our case we are only considering a zero-temperature bath in the case of a damping 
motion caused by friction 



p r (t) = —iQ a^a,p r + ^ (2ap r af — a}ap r — p 



Jr a)a 



The probability of finding n quanta follows a Pauli master equation: 
d 

— (n\ p r \n) = 7 [(n + 1) {n + 1| p r \n + 1) - n (n| p r \n)] . 



(58) 



It is easy to see that p r of Eq. (|50l) , when the background is neglected, is solution of the 
Pauli equation (|58D.H 

The results listed above are well known, but shown that our formalism works as the 
usual well established theories on the subject. 



VI. SEMIGROUPS, WIGNER TIME-INVERSION AND IRREVERSIBILITY 



One of the features of the system we are studying is its irreversible evolution; the system 
reaches the equilibrium at the far future, and of course the inverse evolution is not possible 
anymore. These properties were already found in Refs. |5,21| and here they are briefly 
reviewed. 

The presence of a time-asymmetric behavior can be shown in two different ways: As the 
splitting of the usual evolution group in two semigroups or as the impossibility to make a 
time inversion. We consider that the second one is the most eloquent. 



10 Moreover p r {i) is possitive definites since this condition is equivalent to pupoo > \pio\ 2 (pij 
{i \ p r which is obviously satisfied by p r (t) provided it is satisfied by p r (0). 



22 



A. Semigroups 



Physical states can be chosen to be in test space $ + , as we have done (or, with a simple 
and physically irrelevant change of convention in space $-), while corresponding generalized 
eigenvectors are in its dual (or $*). The proof is simple. From the Paley- Wiener theorem 
H the following lemma can be deduced: 

If f(uj) e Hi, e iujt f(u) e Hi only for positive times. Similarly if f(u) e H 2 , e iujt f(u) G 
H^_ only for negative times. 

The asymmetry in Hardy spaces can be immediately seen from the Fourier transform 
representation. From the Paley- Wiener theorem it is known that Hardy class functions from 
above can be represented as 



/2tt Jo 

where 



f(s) = ^=T duo e-^/M- 

V 27T J— co 



The Fourier transform of Hardy class functions H+ is in the space of square integrable 
functions with support on the positive real axis; and as the Fourier transform of e lu}t f(uj) is 
given by 

r{<**f(u>)} = f{8-t), 

i.e. a function with support on [0, oo) is transformed into a function with support on [t, oo). 
For a negative time this last function has no longer support on [0, oo) and therefore e luJt f(u}) 
does not belong to H+. Analogously it can be proved the same property for the Hardy class 
Hi. 

To simplify, we analyze the one-particle case; generalization to n-particle states is 
straightforward. Let 4>(uj) be a function in 9 S D H± such that, as a consequence of previ- 
ous lemma, e luJt <p(u)) will also be in the same space for t > only. As <f)(uj) = (u}\4>) then we 
can write in Dirac notation 
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e iuJt (j)(uj) = e iujt {to\<P) = (cu\e- iH{ h \(P), (59) 

and taking into account Eq. ( p5|) we can state that if \(f>) E $+ then e~ lH< t \<f>) E $+ only 
for t > 0. In the same way, if \<j>) E $_ then e~ tH(+)t \(p) E for t < 0. Then if we postulate 
that $ + is the space of physical states and that the physical evolution brings physical states 
into physical states, it turns out that this evolution only take place in the period t > 0, so 
irreversibility naturally appears 0. 

B. The Wigner time-inversion 



The Wigner time-inversion operator acts in a real representation as |38 



Ktp(x) = (p*{x). (62) 

[cf. Eq. (H)]. Then, as the complex conjugate of the functions of H\ are the functions of 
H 2 we have that |21 



(63) 



11 Moreover, states in <3?+ are linear combination of generalized vectors \n,uj\...uj ni ) , and these 
vectors evolve as 

e- iH( ~ h \n, LuT^-LOm) = e - i ^ on+ul+ - +u ^ t e~2 nt \ n , uj^.uo m ) . (60) 

Therefore, with the exception of n = 0, they decay towards increasing t. Then we say that physical 
vectors in $ + decay towards positive time (except the states belonging to n as the vacuum 
state 1 0) | v ) , which does not decay) . Analogously 

and we can say that the unphysical states in decay towards negative times. 
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Then the Wigner operator is not well defined either within $ + or within so those 
states in $ + or $_ are not, in general, t— symmetric. Therefore if we consider that only 
the states of $+ are "physical" or "admissible" the Wigner time inversion transforms these 
states in "unphysical" ones and therefore it turns out to be impossible since unphysical states 
simply do not exist in nature. Then, through this mathematical structure, irreversibility is 
incorporated in our theory. 



Nevertheless ($ + R $_) is not an empty set [fLi ], so the time-inversion operator will be 
well defined there 

K : ($ + n $_) — > ($+ n $_), (64) 

and these states will describe reversible processes, which will be t— symmetric. 

All these things lead us to the postulate of the introduction: "Physical states are in $ + 
( or $_)." In fact, this postulate provides a mathematical structure to deal with irreversible 
process, as it was also shown in papers pl,5]-fl|. The choice between $ + and $_ is conven- 



tional and does not lead to physical consequences, but once we choose one of these spaces the 
distinction between past and future becomes substantial. Moreover, if we take into account 
the global structure of the universe, this choice can be motivated from the asymmetry of 
this structure. ( fli~2fl, Ol) 



VII. CONCLUDING REMARKS 



We outline the main results of this work. 

We have diagonalized the full Hamiltonian of our model and extended it in such a way 
that the solution is analytic in the interaction parameter. A rigorous mathematical formalism 
can be introduced in order to deal with unstable quantum systems (see the appendices). 
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Using this formalism, we have obtained a second quantized version of the decay of un- 
stable systems and we have found the corresponding creation and annihilation operators of 
unstable states. 

By means of a simple example, the exact time evolution of the reduced density matrix 
at zero temperature has been studied. We have obtained an exponential decay approach 
of P(t) to the asymptotic value P(t — > oo) = which is expected when the particle is in 
thermal equilibrium with a zero temperature bath. For short times we have found a quadratic 
behavior for the decay probability P(t) (Zeno effect). This short time deviation from the 
exponential decay law was recently measured by first time (see Ref. |[39|| ). Other deviations 



from the exponential decay law, in this case for long times, naturally arise in our framework: 
Khalfin effect, which unfortunately are in practice far enough of any observational time scale. 
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IX. APPENDICES 

APPENDIX A: CREATION AND ANNIHILATION OPERATORS 

Consider the annihilation and creation "unsmeared operators," and 6^, respectively, 
that we have used in our calculations. Commonly they are introduced in the mathematical 



framework of quantum field theory by virtue of expressions like [40 



b{<P)= U (uo) b w duo (Al) 



or 
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fet (0) = / 0* ( W ) 6 t ^ 



(A2) 



where 6 (0) and &t (0) are the (well-defined or "smeared" ) annihilation and creation operators 
of the one-particle state G Ti (being H a Hilbert space) and their action is interpreted as 
the annihilation, respectively the creation, of a spectrum localized quanta, represented by a 
Dirac delta centered in the real value u. In this context we find that 



CXJ CXJ 

b j (0) : sym (ft® n ) -> sym (ft® n 



n=0 



n=0 



&t (0)$ = &t(0) (0 o ,0 1; ...,0 r , 



0,& (0)00,^1(0) 01, •••,^(0)0 



ni • • • i 



(A3) 



where 



^ (0) 0n = 6J, (0) [sym (0« (g, 0( 2 ) ® • • • <g> # 
= sym (0^ 1} <g> 0^ 2) ® • ■ ■ ® 0^ n) ® 0) 



(A4) 



and 6(0) = [^(0)] . 

Of course, in the framework of the Hilbert space foundation of quantum mechanics 
"definitions" (|A1|) and (|A2|) are strictly formal. 

It is evident that these equations are analogous to the formal definition of the Dirac delta 

(x) = J <p {oj) 8(x — to) du>, (A5) 

so, as in the case of the Dirac delta, the rigorous meaning of the "unsmeared operators" 
and blj must be found in the rigged Hilbert space formulation of quantum mechanics. 

One way to do this is considering the explicit definitions of b w and b^ as distribution 
valued operators. Being S C 7i C S x a rigged Hilbert space, we have that 



(JXJ (JXJ 

b u : sym (S® n ) -> sym («S® n ) 



n=0 



n=0 



= b w (0o, 01, . . . , 0„, . . .) = (&lu,01, &2OJ02, • • • , &n^0n, • ■ ■) 



(A6) 



where 



27 



&*A = b nu} [sym (<$ <g> 4>f ® • • • <g> ^ n >) 



(A7) 



and 6 na; acts as in its definition given in Sec. II. Also 

oo r oo 

bl : sym (S® n ) -> sym (<S^) 



n=0 



Ln=0 



= b t (00, 01, • • • , 0„, • • •) = (0, 6 Ow O , &L01, • • • , &L0n, • • •) 



(A8) 



where 



= S?/m g, 0(2) . . . ^ ^{n) g, 



A. 



(A9) 



Observe that, with these definitions, only the normal product is well defined and that 
the "T" symbol is merely a convenient notation, that is, there is not a rigorous Hermitian 
conjugation involved. But if we want to define the canonical commutation relations we will 
be in troubles, because the product of functionals is not uniquely defined. So this is not the 
right way either. 

Fortunately, there is another way to define the annihilation and creation "operators," 
and b^. This point of view is more abstract than the previous one. Remember that the Dirac 
delta can be considered as a tempered distribution, i.e. a continuous linear functional on 
the one-particle regular states space S. In an analogous way the annihilation and creation 
"operators," b u and bj , and all the respective "products" can be considered as continuous 
linear functionals on the Canonical Commutation Relations Algebra = CCR(S). 

We summarize some properties which characterize this kind of algebras JIT] . Remember 
that, being S a nuclear metrizable space, there exists a non- decreasing basis {p a } aeI of 
continuous seminorms such that each seminorm is Hilbertian. Let us denote by 7i a the 
Hilbert space which is the completion of the quotient space S/Ker (p a ) with respect to the 
quotient norm p a = p a /Ker (p a ) , i.e. the space of equivalence classes defined by 

<fi a = {x e S : p a - x) = 0} , 
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where 



p a : S/Ker (p a ) -> R+, 

Pa {<f>a) = Pa (0) • 

The *-algebra CCR(S) is defined as the Hausdorff projective limit (|2| of the collection 



{CCR (H a )} aeI , where CCR (7~t a ) is the C*-algebra |43| generated by the family of operators 
{b (0) : cf) G H a } , with respect to the mappings that inject each CCR(H a ) into CCRil-ip) 
if a > f3, where the order in / is the induced by the ordering of the basis of seminorms. 
We can characterize the CCR (S) as follows. Since every projective limit of a collection of 
C*-algebras is a b*-algebra, i.e. a complete symmetric *-algebra whose topology is defined 
by a basis of continuous submultiplicative seminorms, CCR (S) is also a b*-algebra (see Ref. 
||44j| ) . Moreover, we have that CCR(S) is the strict inductive limit |E| of the collection 
{©j=o s y m | , so CCR (S) is a nuclear strict inductive limit of a collection of Frechet 
spaces or £jF*-algebra |45||. With this we have that the algebra is complete, barreled, and 



nuclear. 

Finally we have that CCR(S) C CCR(H) C [CCR(S)] X , which represents a gener- 
alized Gel'fand triplet. So, viewing the relations (|A1|) and ( |A^ ) as generalized expansions 
in the sense of the well known Maurin's theorem, one can identify and b^ as continuous 
linear functionals on the algebra CCR(S) , and we can say that they belong to [CCR (S)] x . 

APPENDIX B: RIGGED EXTENSION 

In this appendix we find a state space $ with the required properties to implement 
our formalism of unstable states. In order to adequately define the vectors obtained in 
Sec. Ill we must restrict the Hilbert space, which is the basic mathematical structure of 
ordinary quantum mechanics. Recall that Dirac's bras are defined as linear functionals on 
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kets of These functional belong to H* , the topological dual of H. But in this case 
1~C is isomorphic to H, then one works indistinguishably with kets and bras. However, 
if we restrict the topology in order to take a dense subset $ of 7i, we break the one-to- 
one correspondence between elements 4> of $ and continuous linear or antilinear functionals 
F over them. We will call $' the dual space of linear functions and $ x the dual space of 
antilinear functional. We usually use the latter one. It leads to a triplet structure symbolized 
as $ C T~i C $ x , where, to assure the convergence in the norm which defines the topology 
of $, we must require that (<p\F) would be finite @,[y],f|(|. This space $ will be the space 
of "regular" states $ + , as we have explained above. Changing the our convention it can be 
$_ 

In our case, a necessary condition for \<p) e <3> is that, the following expression, a gener- 
alization of Eq. (p5[), would have a rigorous meaning, 



\n,v) = U 



;(-)* 



[a'{z )\2 

POO 



0,V 



Jo [Zo ~ 0Jl\ 



poo poo 

+X n / . . . / dui . . . duj r , 
Jo Jo 



\n-l,ui) + ... 



(Bi) 



Ou;i...u; n ) 



'[z -wi]+ [z -uj n ] + 
The last term of the second member of Eq. (|B1|) must be well defined, so the function 
(<f)\0,Ui...ou n ) = <pQ(ui...u n ) has to have an analytic continuation in each variable Ui (0 < 
% < n) to a region which includes the singularity z Q , so that the integral defines an analytic 
n-dimensional function evaluated at z . 

The simplest choice for (</)\0, 0J\...u v ) which does not depend on the localization of zq is 
that (<j)\0,ui...u n ) would be a Hardy function from below H 2 _ JL4J] for each variable. It is 
equivalent to 



|0,wi...w n ) e 9(SnH_ 



2\®n 



(B2) 



H is the Hilbert space of the states we are considering. It can be the whole space of states or 
some subspace with precise physical properties, as the incoming or outgoing spaces. 
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This generalizes the criterium previously used for the one-particle sector [|5 , 21), 11]. From 
this criterium it can be proved that all the mathematical expressions above are well defined 
(see Appendix C). Then in order to (i, U\...u) m \(f>) {i + m = n, n G M) be well defined, it 
must belong to the following function space 



{l,U) 1 ...U) r , 



m=0 



(B3) 



where S is the Schwartz space |1^J40| , and 9 is the Heaviside step function, which gives the 
restriction to the positive real axis. 

If we do the same in order to define ((f)\rT/v), we find another realization space 



(i,u x ..M m \4>) e($e[snH 2 

m=0 



with (i + m = n, n G Af) . Therefore we define the following spaces 



(B4) 



m=0 



(B5) 



The generalized eigenvectors belong to the dual spaces since they are antilinear 
functionals || on spaces flB3| ) and (|B4[) , 



\i,ui...u m ) G $ 



X 

+ ' 



(B6) 



i,uj l ...u m ) G $„. 



These generalized eigenvectors fulfill the following relations (see 



(i,u 1 ...uj rn \i , J uj' 1 ...u' m ^ 



(B7) 



i,ui...u m \i' 1 u' 1 ...u' m , 



Eqs. 



i, ui...u m \i', uj[...uj' ml ) = 5n '^' EE 5 (< - w n) - 8 (< m " ^ 

\" L -) aeGpreGp 

are valid except for i — i' — 0, and in this case we have 



(B8) 
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0,ui...u n \0,u' 1 ...u' n 



(0,Wi..M) n \Q,u/ 1 ..M/ n ,) 



n i 



\)2 



(B9) 



inl 



\)2 



reG v 



Gp is the group of permutations. Eqs. ( [B7D say that the norm of generalized eigenvectors is 
zero (except i = 0), which is a necessary fact to conserve energy [pi J)|. It is not contradictory 
to have null-norm vectors because these are generalized vectors which are not in the usual 



Hilbert space and have an underlying indefinite metric structure [^TJ . If we define the spaces 
7i + and as the spaces $_ and $ + where the condition about S is not required, we have 



Ti 



± 



(f) I {n,LQ x ...U) m ...\ 



oo 
m=0 



then we arrive to the triplets under Eq. @ and in Eq. (pTOQ .Pl 

Using Eq. (B5| ), we find a double structure of rigged Hilbert spaces for our model, 



$+ c n + c $i 



+• 



(BIO) 



APPENDIX C: THE DOUBLE INTEGRAL THEOREM 



In Ref. |TI] it was demonstrated that if we want that the integral 

g(u)((/>\u) 



L 



duj- 



/R+ Zq — UJ 

would be well defined, it is sufficient that 



(CI) 



\u) G 6(SnH 2 _). 



(C2) 



13 In this way $± is dense in TL± which is the outgoing (incoming) space [46]. The TL cited in 
paper || is actually the outgoing space TL + . 
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In the two variables case it is the integral 

R+ Zq — UJ JR+ 



duo 



,g(uj')((f)\uj,uj') 



Z - UJ 

the one that must be well defined. In this case we prove the following theorem. 
Theorem. 
If 

<f>(u, uj') = (<f)\uj, uj') e 9{S n #!f 2 , 

then integral ( |C3|) is well defined. 
Proof. 

If condition (Pf) is fulfilled, as S is a Frechet space[^]we have ( [Q, page 459) 



(C3) 



(C4) 



(w,wO = X;A i 0iH^(u; , ) J 



i=0 



, 01, 02 



where £^ M < MiH, e n # 2 ) (» = 1,2, . 

the r.h.s. of Eq. (|C5|) is absolutely convergent, namely the series 



(C5) 

when i — > oo, and 



i=0 



is convergent for any continuous seminorm p over 9(S n ifi)® 2 . Let us now define the 
seminorm p 2Q as 



do;- 



da; 



,|<?(u/)||(#,v>| 



(C6) 



'R+ |2 — UJ\ JR+ \Zq — W' 

where D is the distance from R + to z (D = 7/2). We must demonstrate that p ZQ is a 
continuous seminorm. We use the Holder inequality [|4"7| 



z — oj Zq — uo 
for any ^ G C_ (the lower half-plane) .0 Since 



(z - w)(z - W) 



(C7) 



A Frechet space is a metrizable complete space. 
'We remember that 
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Oo - u){zo - to') 
Equation (CI) reads 



sup 



(z - lo){zq - u') 



;lu,uj'e R + > = D 



-2 



(C8) 



(C9) 



for any z G C . Then p zo (4>) is not only a continuous seminorm but also a continuous norm 
over 6(SnH 2 __) m . Then 

g(u)\ f AgWM{u)M(fS)\ 



bo - w'l 



£> 2 V / dw / dw' 



z — oj z — u' 



< oo. 



(CIO) 



So, calling 



2 — ^ z — c</ 

from the corollary of the Lebesgue theorem ( [HH], page 33) we know that, if the series 



(Cll) 



i=0 



converge a.e. in R + x R + , then the series, considered as a function of (u,u'), belongs to L\ 
and 



„ « OO CO „ „ 

/ du,/ du/£./f (u^')=E/ <W du>'f*\LO,u' 



Then, going back to Eq. ( C3|) we have 



i - / |/M|dw, 



(C12) 



sup{|/(u;)|;a; G R+}, 



and also that 5 GLi , i.e. any Schwartz function is integrable. 
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R+ 2 — Co" JR+ Zq — CO JR+ Z — UJ JR+ 



oo 



i=0 



= f A( / b MM / (C13) 

The l.h.s. of Eq. ( |C13| ) is well defined since it is an integral of a L\ function. Moreover it 
is the sum of products of two well defined integrals, like the one of Eq. (JC1|), since from 
hypothesis 4>\(uj), 4>\ (u/) G 9(S PI i/ 2 ). Thus, the proof is complete. 

Of course, this theorem can be generalized from the case of two factors, to the case of n 
factors and, taking into account Eq. (|27|) , it proves that Eq. ( [B"T1) is well defined if condition 



(gH) is fulfilled. 



APPENDIX D: FIGURE CAPTIONS 

Fig. 1: Deformation of the contour of integration taking into account the presence of 
the complex pole z . 

Fig. 2: Behavior of the decay probability showing the Zeno, exponential, and Khalfin phases. 
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